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1. Introduction 

The FIRE marine stratocumulus IFO 
held in July, 1987 produced a data set that 
is far more comprehensive than data sets 
from previous stratocumulus experiments. 
One exciting new development was the use 
of the 10.6/xm lidar system for cloud-top 
mapping that was available on the NCAR 
Electra. This system provided a unique look 
at the small scales of the turbulence in the 
clouds, images of the turbulent structures 
that are quantitatively revealed by condi- 
tional sampling (e.g., Khalsa and Greenhut, 
1987). 

The behavior of these updrafts and 
downdrafts is central to the dynamics of the 
stratocumulus-topped marine boundary 
layer. FIRE’s objectives of understanding 
cloud dynamics and how they affect the 
cloud optical depth ~ which, in turn, is the 
crucial factor in determining the clouds* 
albedo — therefore require the investigation 
of these drafts. This poster discusses initial 
results from a simple model capable of simu- 
lating moist, entraining plumes that are sub- 
ject to water phase changes and radiative 
heating and cooling. The results discussed 
here are limited to plumes that are not 
affected by condensation and evaporation 
but are subject to radiative heating and 
cooling. These results correspond, therefore, 
to the "dry cloud" case discussed by Lilly 
and Schubert (1980). The model’s simpli- 
city limits the realism of the results ~ the 
plumes are assumed not to interact, for 
example ~ but the role of radiative processes 


in influencing the plume dynamics is clear. 
Also revealed is the role the plumes play in 
maintaining the cloud-top inversion. 

This poster abstract presents the model 
equations and methodology used, and 
discusses qualitative results. 

2. Model 

The plume model discussed here is 
based on the work of Telford (1966, 1972) 
and is basically a simplification of the model 
discussed by Chai and Telford (1983). The 
main simplification in the present work con- 
cerns the neglect of downdrafts that are 
explicitly associated with the modeled 
updrafts. In the Telford (1972) model of a 
field of convective plumes, updrafts are sur- 
rounded by downward (compensating) flow; 
here, drafts are completely independent. 

In the complete model, the physical 
processes of the drafts include adiabatic 
expansion, lateral entrainment, buoyant 
acceleration (due to virtual temperature per- 
turbations induced by boundary conditions, 
condensation and radiative heating) and 
parameterized friction. The heating from 
condensation is the most non-intuitive 
aspect of the model equations, and this is 
discussed in detail by Chai (1978) and more 
briefly in Chai and Telford (1983). Here, it 
is ignored. 

Dependent variables in the model are 
the plume radius r, vertical velocity to, rms 
turbulent velocity », total density p, and 
water content p w . Because condensation is 
neglected here, the equation for the conser- 


393 



vation of water is not needed, and therefore 
the four equations represent conservation of 
volume, mass, (vertical) momentum and 
turbulence kinetic energy. Each plume is 
assumed to be in a steady state, and there- 
fore these dependent variables are functions 
of height z. Boundary conditions are 
applied at the surface or at the inversion, for 
updrafts and downdrafts, respectively. The 
plumes are embedded in an environment 
specified by mean temperature profiles T(z) 
which, hydrostatically, determines the pres- 
sure distribution p(z ); pressure perturba- 
tions associated with the plumes are ignored. 

The system of equations for the dependent 
variables is first order and highly non-linear; 
it is solved numerically using a Runge-Kutta 
procedure. 

The equations can be written as 

L (r)-f- - F(r) . (l) 

where 

r(z) 

w(z) 

r « s ,( 2 ) • l 2 > 

Assuming a quiescent (non-turbulent) 
environment and using T(z) and p(z) to cal- 
culate an environmental density profile ~p(z) 
allows the right-hand side of (l) and the 
nonlinear operator L to be written as 
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The environmental vertical velocity w 
can be deduced from the continuity equa- 
tion, provided the area of the affected 
environment is specified. Here, it is assumed 
to be so much larger that the plume that w 
-> 0. The parameter A is an entrainment 
coefficient (~0.08, according to Telford, 
1966), and A is a dissipation coefficient, 
taken to be equal to unity, k is the ratio of 
the gas constant to the specific heat at con- 
stant pressure. The non-adiabatic sources of 
mass are denoted by p. 

In principle, the solution to (1) is given 
simply by 

2 

r(*) = f L -1 Fdz; (5) 

Zq 

since | L | = 2 r 4 w 4 p, the problem is well- 
posed except for plumes of vanishing radius 
or density, or plumes that stop. Accord- 
ingly, given an environmental temperature 
profile and boundary conditions, the system 
(l) can be integrated until the vertical velo- 
city vanishes, at which point, to cite an 
example in physical terms, an updraft either 
disappears or becomes a downdraft. 

In practice, integration of (l) must be 
accomplished numerically because of the 
nonlinearity of F and L. This is simplified 
by defining a density perturbation (which is 
also a temperature perturbation, due to the 
hydrostatic assumption and neglect of pres- 
sure perturbations within the plume) 

6 ee (p/p - l) = (T/T-l), (6) 

where T is the plume temperature. Defining 
an entrainment (inverse) length scale 
2Ai 

A = (1+6); a static stability 

rtu 
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S = (g/c p +dT/dz)/T\ and the parameters 
0 = gS/w 2 and a = -(1 -K.)p~ l dp/dz, where 
g and e p are gravitational acceleration and 
the specific heat of air at constant pressure, 
respectively; and the appropriate heating 
rate Q (this is discussed further below), the 
equations take the form 


ildU 

dz 

dw _ 

dz 

dz 


[-A + /?] w ; 

-A 6 - (5 - Q/w){ 1 + 6) ; 


(7a) 

(7b) 

(7c) 


dU 2 ) 

dz 




rw 


(7d) 


The parameter a acts as an adiabatic expan- 
sion term, and 0 represents a sort of Ber- 
noulli effect, through which accelerating 
plumes become smaller. The heating Q is 
related to the radiative flux divergence by 


Q 


1 ad 

~pc p T dz 


( 8 ) 


In deriving (7c), it has been assumed that 
the radiative heating is horizontally homo- 
geneous, a reasonable first-orcer approxima- 
tion for stratocumulus clouds. In the case of 
isolated cumulus, this is clearly invalid, and 
future work with stratocumulus convection 
will also use differential heating rates. It is 
important to stress that radiative heating 
has a distinct influence on the plume’s per- 
turbation temperature (and therefore its 
dynamics, through the coupled system) even 
though the heating rates are assumed to be 
horizontally homogeneous. The results dis- 
cussed in this poster use radiative fluxes cal- 
culated with the parameterizations 
developed by Hanson and Derr (1987). 
Various assumptions about the (implicit) 
cloud liquid water and solar zenith are used 


to investigate the sensitivity of plumes to 
the heating. 

S. Discussion 

The main purpose of this abstract is to 
present the equations of the plume model 
and briefly to discuss the model’s behavior 
in qualitative terms. Numerical solutions 
will be presented on the poster itself. Even 
without integrating the set of equations (7), 
however, it is possible to infer several 
aspects of the model’s behavior from the 
governing equations. 

Note, first, that, with the exception of 
the term -(S-Q/w) on the right-hand-side 
of (7c), the equations are homogenous and, 
further, take the general form 

= F(f) , (9a) 

where f and F are analogous to T and F in 
(2) and (3). This suggests, in cases of neu- 
tral stability* with no heating, that small 
departures from boundary conditions can by 
deduced analytically with 

f(*) = f 0 exp[-F 0 (2-2 0 )J , (9b) 

where ( ) 0 is a boundary value. For realis- 
tic plumes, | 6\ «1 and | u?| <«, and w and 6 
can take both positive and negative values. 
The behavior of the solutions (9b) can then 
be classified entirely according to the ener- 
getics of the plume in question. Energy- 
producing plumes "" that is, direct circula- 
tions ~ behave identically, be they (warm) 
updrafts or (cold) downdrafts. Similarly, 
energy consuming plumes {indirect circula- 
tions) exhibit the same symmetry. 

The behavior of a warm updraft that 
encounters an abrupt increase in stability is 
also apparent from Eqs. (7). The term 
-5(1 +s5) on the right-hand-side of (7c) will 
immediately cause the temperature pertur- 
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bation to decrease, eventually to become 
negative. When this happens, changes 
sign and the plume decelerates and widens 
[see (7a) and (7b)]. The point at which 
w -> 0 signals the termination of the integra- 
tion. 

Note also that the effect of radiative 
cooling on updrafts will be exacerbated 
tremendously by this behavior within an 
inversion. As the plume slows, the term 
Q/w in (7c) will become disproportionately 
large and will apply positive feedback to the 
the deceleration. The role of radiative cool- 
ing in decelerating updrafts in the presence 
of clouds may thus explain the relatively 
sharper inversion structure observed atop 
the cloud-topped boundary layer when com- 
pared to the un-cloudy boundary layer. In 
the latter case, updrafts are slowed by lack 
of positive buoyancy, and this causes them 
to spread out and stop. In the cloud- topped 
case, the radiative cooling |icts as a braking 
mechanism for the updrafts, stopping them 
nearer the inversion base. The net effect of 
this on many plumes is a thinner inversion 
layer. 

Finally, it can be seen that, in cool 
downdrafts, this behavior works in the oppo- 
site sense, to accelerate the downdraft, but 
with negative feedback. This suggests that, 
near the cloud top, at least, the role of ener- 
getic downdrafts should be much stronger in 
stratocumulus than in the clear boundary 
layer. It also suggests that solar heating will 
affect the different drafts in quite different 
fashions. It is this assymetry of updrafts 
and downdrafts of the same sense of circula- 
tion (i.e., both direct) in the presence of 
radiative heating that causes the dynamics 
of stratocumulus convection to differ so 
dramatically from the clear boundary layer 
heated from below. Further discussion, with 
quantitative results of the cases discussed 
here, will be included on the poster. 
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